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The properties of the spectrum of the overlap Dirac operator and their relation to random matrix theory are 
studied. In particular, the predictions from chiral random matrix theory in topologically non-trivial gauge field 
sectors are tested. 



An important property of massless QCD is 
the spontaneous breaking of chiral symmetry. 
The associated Goldstone pions dominate the 
low-energy, finite-volume scaling behavior of the 
Dirac operator spectrum in the microscopic 
regime, 1/Aqcd « L « l/m w , with L the 
length of the system This behavior can 

be characterized by chiral random matrix the- 
ory (RMT). The RMT description of the low- 
energy, finite-volume scaling behavior is speci- 
fied by symmetry properties of the Dirac oper- 
ator and the topological charge sector being con- 
sidered . The RMT predictions are universal 
in the sense that the symmetry properties, but 
not the form of the potential matters [Q. Fur- 
thermore, the properties can be derived directly 
from the effective, finite-volume partition func- 
tions of QCD of Leutwyler and Smilga, without 
the detour through RMT |§, though RMT nicely 
and succinctly describes and classifies all these 
properties. The topological charge enters the 
RMT prediction via the number of fermionic zero 
modes, related to the topological charge through 
the index theorem. The symmetry properties of 
the Dirac operator fall into three classes, corre- 
sponding to the chiral orthogonal, unitary, and 
symplectic ensembles J3|. Examples are, respec- 
tively, fermions in the fundamental representation 
of gauge group SU(2), fermions in the fundamen- 
tal representation of gauge group SU(iV c ) with 
N c > 3, and fermions in the adjoint representa- 
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tion of gauge group SU(A C ). 

The classification according to the three RMT 
ensembles is connected to the chiral properties of 
the fermions A good non-perturbative regu- 
larization of QCD should therefore retain those 
chiral properties. Until recently such a regular- 
ization was not known. The next best thing were 
staggered fermions, which at least retained a re- 
duced chiral-like symmetry on the lattice. Indeed, 
staggered fermions were used to verify predictions 
of chiral RMT, albeit with two important short- 
comings: (i) staggered fermions in the fundamen- 
tal representation of SU(2) have the symmetry 
properties of the symplectic ensemble, not the or- 
thogonal ensemble as continuum fermions, while 
adjoint staggered fermions belong to the orthog- 
onal ensemble, not the symplectic one. (ii) stag- 
gered fermions do not have exact zero modes at 
finite lattice spacing, even for topologically non- 
trivial gauge field backgrounds, and thus seem to 
probe only the v = predictions of chiral RMT. 

The development of the overlap formalism for 
chiral fermions on the lattice || led to the mass- 
less overlap Dirac operator, a lattice regulariza- 
tion for vector-like gauge theories that retains 
the chiral properties of continuum fermions on 
the lattice Q. In particular, the continuum pre- 
dictions of chiral RMT should apply. Overlap 
fermions have exact zero modes in topologically 
non-trivial gauge field backgrounds [Q, allowing, 
for the first time, verification of the RMT predic- 
tions in v ^ sectors. The nice agreement we 
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shall describe further validates the chiral RMT 
predictions and strengthens the case for the use- 
fulness of the Overlap regularization of massless 
fermions. 

The massless overlap Dirac operator Q is given 

by 

D = ^[l + l5 e(H w (m))} . (1) 

Here, 75-ffu,(— to) is the usual Wilson-Dirac oper- 
ator and e denotes the sign function. The mass m 
has to be chosen to be positive and well above the 
critical mass for Wilson fermions but below the 
mass where the doublers become light. We are 
interested in the low lying eigenvalues of the her- 
mitian operator H = 75 D described in Ref. (?]]. 
We will use the Ritz algorithm Q applied to H 2 
to obtain the lowest few eigenvalues. The numeri- 
cal algorithm involves the action of H on a vector 
and for this purpose one will have to use a rep- 
resentation of e(H w (m)). We used the rational 
approximation discussed in Ref. M] . 



We computed the distribution of the lowest ly- 
ing eigenvalue of the overlap Dirac operator in the 
fundamental representation on pure gauge SU(2) 
configurations with (3 = 1.8 as an example of the 
chiral orthogonal ensemble, on pure gauge SU(3) 
configurations with (3 = 5.1 as an example of the 
chiral unitary ensemble, and in the adjoint rep- 
resentation on pure gauge SU(2) configurations 
with f3 — 2.0 as an example of the chiral sym- 
plcctic ensemble. The lattice size was 4 4 in all 
cases. Chiral RMT predicts that these distribu- 
tions are universal when they are classified ac- 
cording to the three ensembles and according to 
the number of exact zero modes v within each 
ensemble and then considered as functions of the 
rescaled variable z = EyA m i n . Here V is the vol- 
ume and E is the infinite volume value of the chi- 
ral condensate (i/jtp) determined up to an overall 
wave function normalization, which is dependent 
in part on the Wilson-Dirac mass to. RMT gives 
the distribution of the rescaled lowest eigenvalue. 
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A collection of the necessary formulae for the dis- 
tribution of the lowest eigenvalue, P m in{z), can 
be found in ||. 

We compare the RMT predictions with our 
data in Fig. ||. If S is known, the RMT pre- 
dictions for Pmin(z) are parameter free. On the 
rather small systems that we considered here, we 
did not obtain direct estimates of E. Instead, 
we made one-parameter fits of the measured dis- 
tributions, obtained from histograms with jack- 
knife errors, to the RMT predictions, with E the 
free parameter. Our results and some additional 
information are given in Table |. We note the 
consistency of the values for S obtained in the 
v = and v = 1 sectors of each ensemble. Alter- 
natively, we could have used the value of E ob- 
tained in the v = sector, to obtain a parameter 
free prediction for the distribution of the rescaled 
lowest eigenvalue in the v = 1 sector. Obviously, 
the predictions would have come out very well. 

With the fermions in the fundamental repre- 
sentation, we found 81 (for SU(2)), and 147 (for 
SU(3)) configurations with two zero modes and 
1 and 3 with three zero modes. For the orthog- 
onal ensemble, we are not aware of a prediction 
for P m in(z) in the v = 2 sector, while for the uni- 
tary ensemble our data, albeit with very limited 
statistics, agrees reasonably well with the param- 
eter free prediction with E from Table |l|. 

For fermions in the adjoint representation, we 
keep only one of each pair of degenerate eigen- 
values so v — 1 is the sector with two exact zero 
modes. Such configurations cannot be assigned 
an integer topological charge since integer charges 
give rise to zero modes in multiples of four JlO[ , 
and we note there are a significant number of con- 
figurations with two zero modes as seen in Ta- 
ble [j] The good agreement with the RMT pre- 
diction found in this case lends further support 
to the existence of configurations with fractional 
topological charge jl(J ■ 

We have tested the predictions of chiral random 
matrix theory using the overlap Dirac operator 
on pure gauge field ensembles. We find the dis- 
tribution of the lowest eigenvalue in the different 
topological sectors fits well with the predictions of 
chiral RMT, with compatible values for the chiral 
condensate from the different topological sectors. 



Table 1 

The chiral condensate, E, from fits of the distribu- 
tion of the lowest eigenvalue to the RMT predic- 
tions. The third column gives the Wilson-Dirac 
mass parameter used, the fourth the number of 
configurations, N u , in each topological sector. 
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254 


0.2931(45) 



This research was supported by DOE con- 
tracts DE-FG05-85ER250000 and DE-FG05- 
96ER40979. 

REFERENCES 

1. H. Leutwyler and A. Smilga, Phys. Rev. D46 
(1992) 5607. 

2. E. Shuryak and J.J.M. Verbaarschot, Nucl. 
Phys. A560 (1993) 306. 

3. For a recent rev iew, see J.J.M. Verbaarschot, 
hcp-th/9902394 

4. G. Akemann, P.H. Damgaard, U. Magnea and 
S. Nishigaki, Nucl. Phys. B487 (1997) 721; 
M.K. Sener and J.J.M. Verbaarschot, Phys. 
Rev. Lett. 81 (1998) 248. 

5. R. Narayanan and H. Neuberger, Nucl. Phys. 
B443 (1995) 305. 

6. H. Neuberger, Phys. Lett. B417 (1998) 141. 

7. R.G. Edwards, U.M. Heller and R. Naraya- 
nan, Nucl. Phys. B540 (1999) 457. 

8. B. Bunk, K. Jansen, M. Liischer and H. 
Simma, DESY-Report (September 1994); T. 
Kalkreuter and H. Simma, Comput. Phys. 
Commun. 93 (1996) 33. 

9. R.G. Edwards, U.M. Heller, J. Kiskis, and R. 
Narayanan, Phys. Rev. Lett. 82, 4188-4191 
(1999). 

10. R.G. Edwards, U.M. Heller and R. Naraya- 
nan, Phys. Lett. B438 (1998) 96. 



